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£vq , Abstract. We review our current knowledge of the solar chemical composi- 

f^) ■ tion as determined from photospheric absorption lines. In particular we describe 

t-H ' the recent significant revisions of the solar abundances as a result of the appli- 

cation of a time-dependent, 3D hydrodynamical model of the solar atmosphere 
■ instead of ID hydrostatic models. This has decreased the metal content in the 

solar convection zone by almost a factor of two compared with the widely used 
compilation by Anders & Grevesse (1989). While resolving a number of long- 
I ' standings problems, the new 3D-based element abundances also pose serious 

challenges, most notably for helioseismology. 
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' 1- Introduction 

c3 ; 

The chemical composition of the Sun is one of the most important yardsticks 
in astronomy with implications for basically all fields from planetary science 
to the high-redshift Universe. Standard practise is to compare the element 
content in a cosmic object with the corresponding value for the Sun using the 
normal logarithmic abundance scale [X/H] = log (Nx/Nu)^ — log(Nx/Nn) & = 
logex,* — logex,©- Abundance ratios such as [C/Fe] are denned correspondingly. 
Hydrogen is the chosen reference element since it is the most common element in 
the Universe as a whole and that it normally dominates the continuous opacity 
in the formation of the stellar spectra. The solar abundances are, however, 
not only of interest as a reference point in astronomy but are also of profound 
importance for the understanding of our own solar system and the solar interior. 

Most of our knowledge of the overall abundance distribution of the ele- 
ments in the solar system originates in two quite different sources: the solar 
photospheric spectrum and pristine meteorites, in particular the so-called CI 
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chondrites 1 (e.g. Anders & Grevesse 1989). Each has their own advantages and 
disadvantages. The abundances of a very large number of elements can be deter- 
mined using absorption lines in the solar spectrum. However, since the results 
are dependent on having realistic models of the solar atmosphere and the line 
formation process the accuracy can only be expected to be good at the 10% 
(0.04 dex) level at best. Isotopic abundance information is only available in very 
rare instances. Model-independent meteoritic isotopic abundances can be mea- 
sured to exquisite accuracy in the laboratory but only a very small subsample of 
meteorites do not show effects of fractionation. In addition, elements like H, He, 
C, N, O and Ne are all volatile and hence depleted in the meteorites. It is there- 
fore clear that the two methods complement each other well. The present review 
will describe our knowledge of the solar photospheric chemical composition, in 
particular in light of recent significant revisions of these abundances following 
the application of the new generation of three-dimensional (3D) hydro dynamical 
model atmospheres to the problem. 

2. Model Atmospheres, Line Formation and Other Obstacles 

While the chemical composition of a star is encoded in its spectrum, it is not 
immediately extractable without first modelling the stellar atmosphere and the 
line formation process. Since stellar element abundances are thus not observed 
but derived they are by necessity model-dependent. As more realistic modelling 
become available the estimated stellar abundances may therefore change with 
time. Equally important is improved atomic and molecular data, which is nicely 
exemplified by Sneden et al. in these proceedings. Obviously progress in both 
of these areas is vital. Indeed they are also complementary as more realistic 
modelling can uncover shortcomings in the input physics and vice versa. 

Traditionally, standard stellar abundance analyses employ time-independent, 
one-dimensional (ID) model atmospheres in hydrostatic equilibrium. Such mod- 
els basically come in two types: theoretical and semi-empirical. In the former 
the temperature structure is computed assuming flux constancy, i.e. the same 
energy is transported through each atmospheric layer. The energy is assumed 
to be carried solely by radiation or convection. For simplicity, the convective 
energy flux is approximated with the local mixing length theory (Bohm-Vitense 
1958) or some close relative thereof (e.g. Canuto & Mazzitelli 1991). In the 
absence of convection, radiative equilibrium holds. It is crucial to include the 
effects of blanketing from spectral lines when integrating the radiative transfer 
equation over all wavelengths in order to obtain a realistic temperature struc- 
ture. LTE is normally assumed for the computation of the atomic populations 
and the radiative transfer. Some of the most widely-used theoretical ID model 
atmospheres for late-type stars are those by Kurucz (1979), the MARCS consor- 
tium (Gustafsson et al. 1975) and the Phoenix team (Hauschildt et al. 1999), 
which are all continuously improving their models, mainly through the inclusion 
of better and more complete atomic and molecular physics. 



1 Important and complementary information is provided by studies of for example the solar 
corona, solar wind, solar energetic particles, solar flares and lunar soil, which however will not 
be discussed here. 
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In a semi-empirical model atmosphere the temperature is deduced from 
observations through an inversion procedure, normally using spectral lines of 
different line formation depths and/or center-to-limb variations. It therefore 
follows that the atmosphere in such a model does not need to obey flux constancy 
nor is it necessary to estimate the radiative and convective energy transport. 
The danger with such a method is that the line formation calculations may 
be faulty, for example due to the assumption of LTE or the use of erroneous 
line data. This would directly translate to an incorrect temperature structure. 
Hydrostatic equilibrium is assumed also in semi-empirical model atmospheres. 
Likewise, a knowledge of the equation-of-state is necessary to relate temperature, 
gas pressure and electron pressure. Semi-empirical atmosphere modelling has 
only been tried relatively few times for stars other than the Sun (e.g. Allende 
Prieto et al. 2000) but is a relatively common practise in solar analyses. Two 
of the most widely used solar semi-empirical models are the VAL-3C (Vernazza 
et al. 1976), and Holweger-Miiller (1974) models. The former is optimised for 
chromospheric studies while the latter has traditionally been the model of choice 
for solar abundance analyses. 

More recently, a new generation of time-dependent 2D and 3D hydrodynam- 
ical model atmospheres has become available and applied in stellar abundance 
analyses (e.g. Nordlund &: Dravins 1990; Dravins & Nordlund 1990; Steffen et al. 
1995; Atroshchenko & Gadun 1994; Stein & Nordlund 1998; Asplund et al. 1999, 
2000a,b, 2004a,b; Asplund 2000, 2004; Asplund & Garcia Perez 2001; Steffen 
& Holweger 2002; Vogler et al. 2004). By simultaneously solving the standard 
hydrodynamical equations for conservation of mass, momentum and energy to- 
gether with the 3D radiative transfer equation, both the radiative and convective 
energy transport are self-consistently computed. No mixing length theory with 
its free parameters therefore enters the modelling. The effects of line-blanketing 
are accounted for through the use of opacity binning assuming LTE (Nordlund 
1982). State-of-the-art equation-of-state (e.g. Mihalas et al. 1988) and opacities 
(e.g. Kurucz 1993; Cunto et al. 1993) are used. One can easily imagine that 
for late-type stars such as the Sun where the convection zone reaches up to and 
even beyond the optical surface, such 3D hydrodynamical modelling is neces- 
sary to capture the time-dependence and resulting atmospheric inhomogeneities 
of granulation. Being more sophisticated does not, however, necessarily make 
such 3D models more realistic. A great deal of effort lately has therefore been 
invested in testing the predictions of these models against an arsenal of observa- 
tional constraints, in particular for the Sun. So far the comparison has been very 
satisfactory, which gives confidence in the realism of the 3D models and their 
predictions (e.g. Stein & Nordlund 1998). For example, the intrinsic Doppler 
shifts inherent in the hydrodynamical modelling render the classical concepts of 
micro- and macroturbulence obsolete, yet the theoretical 3D line profiles agree 
essentially perfectly with the observed profiles, including their line shifts and 
asymmetries (e.g. Asplund et al. 2000a). 

Besides the atmospheric modelling, it is also necessary to accurately model 
the line formation and the radiative and collisional processes that couple the 
radiation field with the gas. Most analyses are based on the assumption of LTE, 
which simplifies matter immensely by assuming that the atomic populations are 
described fully by the Boltzmann and Saha distributions for excitations and 
ionizations. In reality, for conditions typical of stellar atmospheres, the colli- 
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sional rates do not dominate over the corresponding radiative rates, which is 
a pre-requisite for LTE to hold. Often it is therefore important to account for 
departures from LTE, which can be substantial in terms of derived abundances. 
The numerical techniques and codes to compute the statistical equilibrium in 
ID have been available for decades, indeed even 3D non-LTE calculations are 
now feasible (e.g. Asplund et al. 2003). In spite of this, remarkably little work 
has been devoted to this important aspect of abundance analyses even in ID. 
In fact, to date only a dozen or so elements have had their non-LTE line forma- 
tion in solar-type atmospheres studied in some detail with contemporary atomic 
physics. Perhaps this can blamed on a notion within the astronomy community 
that the input physics is very uncertain. This is, however, only partly true, 
since the last decade has seen tremendous progress in this respect. Further- 
more, it is possible to test the sensitivity of the outcome to the various rates, 
which can constrain the non-LTE effects. It is important to realize that LTE is 
more likely than not to be wrong with departures from LTE often being on the 
level of 0.1 dex for the Sun. Depending on the required accuracy in the derived 
abundances, LTE may or may not therefore be an acceptable approximation. 

3. The Solar Photospheric Chemical Composition 

The solar abundance analysis which is described below has been carried out using 
a 3D hydrodynamical model of the solar atmosphere (Asplund et al. 2000a). 
Both atomic and molecular lines of very different temperature and pressure 
sensitivities have been employed whenever possible in the hope that it should 
minimise systematic errors. A special effort has been made to utilise only the 
best available line data. In several incidences detailed non-LTE calculations 
have been carried out or the results thereof been culled from the literature. 
As a result, we are confident that the end-product is the most reliable element 
abundance analysis of the Sun carried out so far. Table ^presents a compilation 
of the most reliable solar and meteoritic abundances in our opinion. 

3.1. Lithium, Beryllium and Boron 

Lithium is depleted in the solar convection zone by about a factor of 160 and as 
a consequence even the Lii resonance line at 670.8 nm is very weak (equivalent 
width W\ ~ 0.18 pm=l. 8mA); the subordinate line at 610.4 nm is hopelessly 
puny for the purpose of deriving a solar Li abundance. In addition, the 670.8 nm 
line is blended with CN and Fei lines, which makes the determination of the 
solar Li abundance rather uncertain. Furthermore, Lii is a minority species, 
which is highly sensitive to both the temperature and the photo-ionizing UV 
radiation field (Kiselman 1997; Uitenbroek 1998). Our recommended solar Li 
abundance stems from the ID LTE analysis by Miiller et al. (1975). To this 
value we add the effects of atmospheric inhomogeneities and 3D non-LTE line 
formation (Asplund et al. 1999; Asplund et al. 2003; Barklem et al. 2003) and 
correct for a revised transition probability to arrive at logeLi = 1.05 ± 0.10. 

Until recently Be was thought to be depleted in the solar photosphere by 
about a factor of two (Chiemlewski et al. 1975; Anders & Grevesse 1989). Bal- 
achandran & Bell (1998), however, revisited the issue by attempting to estimate 
any missing UV opacity and found a Be abundance from the Be n 313.1 nm reso- 
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Tabic 1. Element abundances in the present-day solar photosphere and in 
meteorites (CI chondrites). Indirect solar estimates are marked with [..] 





Elem. 


Photosphere 


Meteorites 




Elem. 


Photosphere 


Meteorites 


1 


H 


12.00 


a a r" i a at 

8.25 ± 0.05 


44 


Ru 


1.84 ± 0.07 


-\ '~7'~T 1 A AO 

1.77 ± 0.08 


2 


He 


[10.93 ± 0.01J 


1.29 


45 


Rh 


1.12 ± 0.12 


1 A T I A AA 

1.07 ± 0.02 


3 


Li 


1 nr i n -i n 

1.05 ± 0.10 


o n f i a a/ 1 

3.25 ± 0.06 


46 


Pd 


1.69 ± 0.04 


1 C <~~? I AAA 

1.67 ± 0.02 


4 


Be 


"1 O O 1 A A A 

1.38 ± 0.09 


-i o o i a An 

1.38 ± 0.08 


47 


A 

Ag 


0.94 ± 0.24 


-1 AA 1 A i \ / ' 

1.20 ± 0.06 


5 


B 


a t a i a a r\ 

2.70 ± 0.20 


A T r* 1 A A A 

2.75 ± 0.04 


48 


Cd 


1 i nil 

1.77 ± 0.11 


-i ry-i | a ao 

1.71 ± 0.03 


6 


c 


o on i /"i nr 

8.39 ± 0.05 


^— ,( r\ i a A/ 1 

7.40 ± 0.06 


49 


In 


1.60 ± 0.20 


A OA 1 A AO 

0.80 ± 0.03 


7 


N 


T TO 1 A f\C 

7.78 ± 0.06 


r> c\t* i a at 

6.25 ± 0.07 


50 


bn 


c~\ r\/~\ i /n oa 

2.00 ± 0.30 


AAA 1 A A ,( 

2.08 ± 0.04 


8 





o r* r» i a a f 

8.66 ± 0.05 


O OA 1 A A A 

8.39 ± 0.02 


51 


Sb 


-| A A 1 A O A 

1.00 ± 0.30 


"1 AO 1 A AT 

1.03 ± 0.07 


9 


1 


4.56 ± 0.30 


< a . > 1 a a r> 

4.43 ± 0.06 


52 


Te 




A -\ A 1 A A , ( 

2.19 ± 0.04 


10 


Ne 


l~T o A i a a /~"1 

[7.84 ± 0.06J 


-1.06 


53 


I 




1.51 ± 0.12 


11 


AT 

JNa 


r* -i 1-7 | A n 

6.17 ± 0.04 


r> a T | a r\ o 

6.27 ± 0.03 


r a 

54 


Xe 


To n^7 i n nnl 

[2.27 ± 0.02 


-1.97 


12 


Mg 


n o i a a a 

7.53 ± 0.09 


T r - o i r\ ah 

7.53 ± 0.03 


55 


Cs 


1 A f~~7 I A A O 

1.07 ± 0.03 


13 


Al 


/"* OT 1 ("I A 

6.37 ± 0.06 


6.43 ± 0.02 


56 


Ba 


A 1 f7 I A A T 

2. 17 ± 0.07 


A "1 I A A O 

2. 16 ± 0.03 


14 


Ci • 

Si 


7.51 ± 0.04 


7.51 ± 0.02 


57 


La 


-i -i o i n n r* 

1.13 ± 0.05 


1 1 r - in n c 

1.15 ± 0.06 


15 


P 


f » ) / < I A. A , | 

5.36 ± 0.04 


r* a r\ i n n 

5.40 ± 0.04 


58 


Ce 


1 r* o i n nn 

1.58 ± 0.09 


-1 f A 1 A A A 

1.58 ± 0.02 


16 


s 


T -i a I a r\f 

7.14 ± 0.05 


r-r -\ r* i n n A 

7.16 ± 0.04 


59 


Pr 


A f-7-l 1 AAA 

0.71 ± 0.08 


A r—r 1 A AO 

0.75 ± 0.03 


17 


CI 


p" v a i a oa 

5.50 ± 0.30 


r* oo i n n/ 1 

5.23 ± 0.06 


60 


AT 1 

Nd 


1.45 ± 0.05 


1 JO I A A O 

1.43 ± 0.03 


18 


A 

Ar 


[6. 18 ± 0.08J 


-0.45 


62 


Ci 

am 


-\ A 1 1 A A r> 

1.01 ± 0.06 


AAA I A A A 

0.92 ± 0.04 


1 n 

19 


K 


r* r\o i n n t 

5.08 ± 0.07 


5.06 ± 0.05 


63 


Eu 


i\ r- a i AA/ 1 

0.52 ± 0.06 


A A A I A A I 

0.49 ± 0.04 


20 


Ca 


/i o -1 1 A A , ( 

6.31 ± 0.04 


r> c\r\ i /^i no 

6.29 ± 0.03 


f* A 

64 


Gd 


i in i n n ^ 

1.12 ± 0.04 


i no i n nn 

1.03 ± 0.02 


21 


Sc 


O A P* 1 A A O 

3.05 ± 0.08 


3.04 ± 0.04 


65 


Tb 


A AA k A i\ { \ 

0.28 ± 0.30 


A AA 1 A AO 

0.28 ± 0.03 


22 


Ti 


A nn i n n c 

4.90 ± 0.06 


4.89 ± 0.03 


66 


Dy 


-| -| A 1 AAA 

1.14 ± 0.08 


-| -1 A 1 A A 1 

1.10 ± 0.04 


23 


V 


J a a i a a a 

4.00 ± 0.02 


3.97 ± 0.03 


67 


Ho 


A V "1 1 A I A 

0.51 ± 0.10 


A A r* 1 A AA 

0.46 ± 0.02 


24 


Cr 


5.64 ± 0.10 


r* <^ o i /~i at 

5.63 ± 0.05 


68 


Er 


a ao i a r\r> 

0.93 ± 0.06 


A AA 1 A AO 

0.92 ± 0.03 


25 


Mn 


r on i n no 

5.39 ± 0.03 


r ^ t i n no 

5.47 ± 0.03 


69 


lm 


n nn i Air 

0.00 ± 0.15 


n no i n nr 1 

0.08 ± 0.06 


26 


Fe 


7.45 ± 0.05 


t—r A r* 1 /^i AO 

7.45 ± 0.03 


70 


Yb 


"1 A A 1 Alp" 

1.08 ± 0.15 


A A "I 1 A A O 

0.91 ± 0.03 


27 


Co 


4 A A | A A C*i 

4.92 ± 0.08 


A in r* \ r\ r\ o 

4.86 ± 0.03 


71 


Lu 


A A r* 1 A ~\ A 

0.06 ± 0.10 


a a r* i a r\r> 

0.06 ± 0.06 


28 


Ni 


r* c\ o i a a ^ 

6.23 ± 0.04 


z 1 1 r\ i n n o 

6.19 ± 0.03 


72 


Hi 


A AO 1 A AO 

0.88 ± 0.08 


A f~7 A 1 A A A 

0.74 ± 0.04 


29 


Cu 


i a -i i r\ r\ A 

4.21 ± 0.04 


4.23 ± 0.06 


73 


Ta 




A "I l~r 1 A A o 

-0.17 ± 0.03 


30 


Zn 


A r* c\ i a no 

4.60 ± 0.03 


4.61 ± 0.04 


74 


W 


1.11 ± 0.15 


A r> A 1 A AO 

0.62 ± 0.03 


31 


Ga 


2.88 ±0.10 


3.07 ±0.06 


75 


Re 




0.23 ± 0.04 


32 


Ge 


3.58 ±0.05 


3.59 ±0.05 


76 


Os 


1.45 ±0.10 


1.34 ±0.03 


33 


As 




2.29 ±0.05 


77 


Ir 


1.38 ±0.05 


1.32 ±0.03 


34 


Se 




3.33 ±0.04 


78 


Pt 




1.64 ±0.03 


35 


Br 




2.56 ±0.09 


79 


Au 


1.01 ±0.15 


0.80 ± 0.06 


36 


Kr 


[3.28 ± 0.08] 


-2.27 


80 


Hg 




1.13 ±0.18 


37 


Rb 


2.60 ±0.15 


2.33 ±0.06 


81 


Tl 


0.90 ±0.20 


0.78 ± 0.04 


38 


Sr 


2.92 ±0.05 


2.88 ± 0.04 


82 


Pb 


2.00 ±0.06 


2.02 ±0.04 


39 


Y 


2.21 ±0.02 


2.17 ±0.04 


83 


Bi 




0.65 ± 0.03 


40 


Zr 


2.59 ±0.04 


2.57 ±0.02 


90 


Th 




0.06 ± 0.04 


41 


Nb 


1.42 ±0.06 


1.39 ±0.03 


92 


U 


<-0.47 


-0.52 ± 0.04 


42 


Mo 


1.92 ±0.05 


1.96 ±0.04 
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nance doublet in agreement with the meteoritic value. There is a long-standing 
debate whether or not there is a need for additional continuous and/or line 
opacity in the UV spectral region over those already included in the model cal- 
culations in order to properly model the solar flux distribution. The novel feature 
with the analysis of Balachandran &: Bell was their estimation of the amount 
of missing UV opacity by requiring that the OH A-X electronic lines around 
313 nm should yield the same solar oxygen abundance as the OH vibration lines 
in the IR. In order to achieve this, they estimated that about 60% extra opacity 
was required, which they attributed largely to photo-ionization of Fe I (see also 
Bell et al. 2001). A follow-up study using a 3D solar model atmosphere instead 
of the MARCS (Gustafsson et al. 1975) and Holweger-Miiller (1974) ID models 
employed by Balachandran & Bell confirmed that Be is undepleted in the solar 
convection zone: logeee = 1.38 ±0.09 (Asplund 2004). It should be emphasized, 
however, that this result hinges completely on the correctness of the utilised 
method of estimating the missing opacity. In particular, departures from LTE 
for the OH A-X electronic lines could skew the results. No sufficiently detailed 
non-LTE calculation exists for these lines to our knowledge, which obviously is 
important in order to settle the issue. Likewise, improved determinations of the 
Fe I bound-free cross-sections would be very valuable. 

The problem with possible missing UV opacity may also affect the deter- 
mination of the solar B abundance, as the only available reliable indicator is the 
Bi 249.67 nm resonance line. Cunha & Smith (1999) presented a re-analysis of 
the Bi line, including a careful study of possible missing opacity sources such 
as photo-ionization of Mg I. Whilst the uncertainty in their result was relatively 
large, they found no evidence for a solar B depletion, which of course would be 
highly surprising if indeed Be is undepleted. We have estimated the 3D LTE 
effect (i.e. 3D LTE - ID LTE) on the resonance line to be 0.04 dex in compari- 
son with a solar marcs or Kurucz model and —0.04 dex for the Holweger-Miiller 
model. When also including departures from LTE as predicted by Kiselman 
&; Carlsson (1996), we obtain logee = 2.70 ± 0.20. The large uncertainty is 
dominated by the problem in analysing this very crowded spectral region. 

3.2. Carbon 

The solar carbon abundance can be determined from a wide range of indica- 
tors, both atomic and molecular transitions. We have recently performed a 
detailed study of the best abundance diagnostics in our opinion: the forbidden 
[Cl] 872.7 nm line, high excitation Ci lines, CH vibration-rotation IR lines, CH 
A-X electronic lines, C2 Swan electronic lines and CO vibration lines (Allende 
Prieto et al. 2002; Asplund et al. 2004b; Scott et al., in preparation). The 
line formation was computed based on a 3D hydro dynamical model of the solar 
atmosphere (Asplund et al. 2000a) and departures from LTE were accounted 
for in the case of the C 1 lines. In sharp contrast with analyses using ID model 
atmospheres, very gratifying agreement is found between all abundance indica- 
tors when employing the 3D model, as clear from Table |21 Given the extreme 
temperature sensitivity of molecule formation, it is crucial to include the effects 
of atmospheric inhomogeneities for the molecular lines in order not to estimate 
too high abundances, which is the case with ID models. Our new solar C abun- 
dance of logec = 8.39 ± 0.05 is 0.17 dex lower than the recommended value in 
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Tabic 2. The derived solar C, N and O abundances as implied by a variety 
of different atomic and molecular indicators using a 3D hydrodynamical model 
of the solar atmosphere (Asplund et al. 2004a, b; Asplund ct al., in prepara- 
tion; Scott et al., in preparation). In addition, the corresponding results for 
two ID model atmospheres, a theoretical MARCS (Asplund et al. 1997) and 
the semi-empirical model of Holweger-Miiller (1974), are given for compari- 
son purposes. Note the excellent agreement between atomic and molecular 
diagnostics with the 3D model and the large spread when using ID model 
atmospheres. The quoted uncertainty is only the line-to-line scatter; no er- 
ror is given for [Cl] as the result is based only on the 872.7 nm line. The O 
abundance used for the CO calculations was 8.66 in the case of 3D and 8.85 
for the two ID models 



lines 






logec,N,o 








3D 


HM 


MARCS 


[Cl] 




8.39 


8.45 


8.40 


Ci 




8.36 ± 0.03 


8.39 ± 0.03 


8.35 ± 0.03 


CH Av = 


1 


8.38 ± 0.04 


8.53 ±0.04 


8.42 ± 0.04 


C2 Swan 




8.44 ± 0.03 


8.53 ±0.03 


8.46 ± 0.03 


CH A-X 




8.45 ±0.04 


8.59 ±0.04 


8.44 ± 0.04 


CO Av = 


1 


8.41 ±0.02 


8.62 ± 0.02 


8.57 ±0.02 


CO Av = 


2 


8.38 ±0.02 


8.70 ± 0.03 


8.59 ±0.03 


Ni 




7.85 ± 0.08 


7.97 ±0.08 


7.94 ± 0.08 


NH Av = 


1 


7.73 ± 0.05 


7.95 ± 0.05 


7.82 ± 0.05 


[Oi] 




8.68 ±0.01 


8.76 ±0.02 


8.72 ±0.01 


Oi 




8.64 ± 0.02 


8.64 ± 0.08 


8.72 ± 0.03 


OH At; = 





8.65 ±0.02 


8.82 ±0.01 


8.83 ± 0.03 


OH At; = 


1 


8.61 ±0.03 


8.87 ±0.03 


8.74 ± 0.03 



the widely used compilation of Anders & Grevesse (1989). The excellent agree- 
ment between transitions of very different formation depths and temperature 
and pressure sensitivities is a very strong argument in favour of this low abun- 
dance as well as for the realism of the 3D model. In particular we note with 
satisfaction that consistent results are now finally provided also by the CO lines, 
which have previously caused a great deal of trouble when analysed within the 
framework of ID model atmospheres (Grevesse et al. 1995; Ayres 2002 and 
references therein). 

A 3D analysis of weak 12 C 16 and 13 C 16 lines yields a photospheric 
12 C/ 13 C ratio of 95 ± 5 (Scott et al., in preparation), which is in excellent 
agreement with telluric measurements (Rosman & Taylor 1998). 

3.3. Nitrogen 

We have recently re-determined the solar N abundance on the basis of the same 
3D hydrodynamical model atmosphere as previously used for the corresponding 
C and O analyses (Asplund et al, in preparation). In the case of N there are 
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fewer abundance indicators compared with its two neighboring elements. We 
have utilised Ni and NH vibration-rotation lines while relegating the CN red 
band to a supporting role. The Ni list consists of 21 weak (W\ = 0.1 — 1pm), 
high-excitation (xexc = 10.3 — 12.4 eV) lines, which implies a 3D LTE abundance 
of log 6n = 7.88±0.08. No 3D non-LTE study for N is as yet available but the ap- 
plication of the computed ID non-LTE abundance corrections by Rentzsch-Holm 
(1996) lowers this abundance to logeN = 7.85 db 0.08. It should be noted, how- 
ever, that Rentzsch-Holm (1996) employs efficient inelastic H collisions, which 
help thermalise the level populations. Given that the available laboratory and 
quantum mechanical calculations tend to suggest that the classical recipe used 
by her significantly over-estimates the H collision cross-sections (e.g. Belyaev 
et al. 1999; Barklem et al. 2003), the Ni based abundance presented here may 
require a small downward revision. 

We have also used a sample of 22 weak NH vibration-rotation lines in the 
IR (located around 3/xm), which have previously been advocated as primary 
abundance indicators (Grevesse et al. 1990). As with all molecular lines, the 
NH lines are very temperature sensitive and as consequence there is a large 
difference between the results of ID and 3D analyses. The NH lines indicate an 
abundance of log eN = 7.73 ±0.05 with our 3D solar model, which is significantly 
lower than with the standard Holweger-Miiller ID model. 

3.4. Oxygen 

Allende Prieto et al. (2001) argued for a significantly lower solar O abundance 
than the then widely accepted value, largely due to the presence of a previously 
unaccounted for blend in the forbidden [O i] 630.0 nm line. We have subsequently 
performed a much more detailed study, which also included high-excitation per- 
mitted O I lines and molecular transitions of OH, both vibration-rotation and 
pure rotation lines (Asplund et al. 2004a). As the Oi lines are known to be 
susceptible to non-LTE effects, full 3D non-LTE line formation calculations were 
performed for these lines. It has long been clear that the OH lines imply higher 
abundance by about 0.2 dex than the Oi lines when including departures from 
LTE, regardless of which ID model atmosphere has been used. Finally this prob- 
lem has been resolved with the application of 3D hydro dynamical models. The 
agreement between the various atomic and molecular lines is excellent, in spite 
of their very different atmospheric sensitivities. The almost perfect match be- 
tween predicted and observed line profiles and center-to-limb variations provides 
further evidence that the 3D model is highly realistic. Furthermore, Melendez 
(2004) finds a similarly low O abundance from the extremely weak first overtone 
OH vibration lines, which have previously not been used. 

A 3D analysis of isotopic CO lines yields a photospheric 16 0/ 18 ratio of 
530 ± 40 (Scott et al., in preparation), which is in good agreement with the 
telluric data (Rosman &: Taylor 1998). 

3.5. Neon and Argon 

The photospheric abundances of Ne and Ar can not be determined directly due 
to the lack of suitable spectral lines. Their recommended values in Table ^ 
have been estimated from measured abundance ratios in the solar corona and 
solar energetic particles in conjunction with the photospheric abundance for 
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the reference element, which in this case is O. The Ne and Ar abundances are 
therefore directly affected by the revised solar O abundance. This explains most 
of the difference compared with for example Anders & Grevesse (1989) with a 
secondary effect coming from refined coronal abundance ratio determinations 
(Reames 1999). 

3.6. Sodium to Calcium 

In connection with the present conference, we have performed a preliminary 3D 
LTE abundance analysis for the remaining elements up to calcium. We have 
been heavily influenced by the study of Lambert & Luck (1978) who performed 
a similar analysis using the Holweger-Miiller (1974) semi-empirical ID model 
atmosphere. In particular our choice of lines has been largely moulded by their 
recommendations, with the added constraint that we have so far restricted the 
comparison to only weak lines. Line profile fitting has been used in all cases to 
estimate the abundance instead of equivalent widths. We have ransacked the 
recent literature for the most up-to-date transition probabilities. As appropriate, 
we take into account departures from LTE, although computed for ID models 
since full 3D non-LTE calculations for these elements are well beyond the scope 
of the present study. While this is inconsistent in principle it is expected to 
give at least a good first estimate of the 3D non-LTE effects (Asplund et al. 
2004a). Space permits only a very brief summary here with the full details being 
presented elsewhere (Asplund et al., in preparation). The preliminary results 
are listed in Table ^ The new 3D-based abundance are in all cases slightly 
lower than those recommended by Anders &z Grevesse (1989) and Grevesse & 
Sauval (1998), the difference being typically 0.05 — 0.10 dex. The impact of the 
3D model atmosphere is thus smaller here than for C, N and O, mainly since 
these abundances are based on atomic transitions rather than very temperature 
sensitive molecular lines. 

The Na abundance is based on six weak Nai lines. The computed non- 
LTE abundance corrections are small (< 0.05 dex) but slightly lower the mean 
abundance to logeNa = 6.17 ± 0.04, which is 0.1 dex lower than the meteoritic 
value. As the ^/-values should be reliable, an explanation for this minor dis- 
crepancy has not yet been identified. The transition probabilities for Mgi are 
notoriously uncertain, forcing the Mg abundance to be estimated from four Mg n 
lines: logeMg = 7.53 ±0.09. Non-LTE effects should be negligible for Mgn. The 
Mgi value is very similar but with unacceptably large uncertainties. Many All 
lines are detected in the solar spectrum but the majority give quite discrepant 
results, mainly due to cancellation and configuration interaction effects in the 
calculations of the ^/-values. We consider only four Al I lines to be of sufficiently 
high quality for this purpose, which suggests logeAi = 6.37 ± 0.06. Based on 
the published departure coefficients by Baumiiller & Gehren (1996), possible 
non-LTE effects should be minor for these relatively high-excitation lines. The 
Si abundance is that of Asplund (2000). 

Our recommended P abundance stems from an analysis of five P I lines in 
the near-IR with equivalent widths W\ < 2.5 pm. With the ^/-values of Berzinsh 
et al. (1997) we find logep = 5.36 ± 0.04. From the line list used by Lambert 
&; Luck (1978), we removed several Si lines with profiles that indicate blends. 
The pruned list of ten lines indicate loges = 7.14 ± 0.05. The new S abundance 
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is now in excellent agreement with the meteoritic value. Chlorine, like fluorine, 
has not been possible to analyse with our 3D model atmosphere of quiet Sun as 
the relevant HC1 lines are only present in sunspot spectra. The highly uncertain 
CI abundance dates back to Hall & Noyes (1972). As we have restricted the 
analysis to weak and intermediate strong lines, we have not yet used the very 
strong Ki 769.8 nm resonance line and instead base our K abundance on five 
weak K I transitions at 404.4, 580.1, 1176.9, 1243.2 and 1252.2 nm. Together with 
the non-LTE abundance corrections of Ivanova & Shimanskii (2000), these lines 
imply log ex = 5.08 ± 0.07. Finally, for Ca a nice selection of both permitted 
and forbidden Cai transitions as well as Call lines are available in the solar 
spectrum. The agreement between the two ionization stages is excellent, and we 
simply take the average of the two: logeca = 6.31 ± 0.04. 



3.7. Other Elements 

Of the remaining elements, only Fe has been the topic of a 3D abundance analysis 
(Asplund et al. 2000b), whose value we adopt here: logeFe = 7.45±0.05. Almost 
identical abundances are suggested by Fe I and Fe n lines. We note that this value 
has not been corrected for departures from LTE. Shchukina & Trujillo Bueno 
(2001) found that the Fei based abundance should be adjusted upwards by 
about 0.05 dex based on 1.5D (i.e. only including vertical radiative transfer and 
therefore treating each atmospheric column as a separate ID model atmosphere) 
non-LTE calculations using one snapshot from the 3D solar model used here. 
Their calculations did not include any inelastic H collisions, which would tend 
to decrease the non-LTE effects towards restoring the LTE results. The jury 
is still out whether or not these collisions are efficient. Korn et al. (2003) 
present convincing arguments that they are in the case of Fe, and we therefore 
tentatively recommend the 3D LTE value. Clearly, detailed quantum mechanical 
calculations for Fe+H collisions are very important but unfortunately the task 
is extremely challenging. 

While no 3D studies have been carried out for the remaining elements, the 
recommended values in Table^nevertheless differ in several incidences compared 
with previous compilations. This is largely the result of the very important 
work done by a small group of dedicated atomic physicists in measuring and 
calculating improved transition probabilities. Much recent attention has been 
devoted to the heavy elements formed by neutron capture, like La (Lawler et al. 
2001a), Nd (Den Hartog et al. 2003), Eu (Lawler et al. 2001c), Tb (Lawler et al. 
2001b), Ho (Lawler et al. 2004), Pt (Den Hartog et al. 2004) and Pb (Biemont 
et al. 2000). Some beautiful examples of this are shown in Sneden et al. (these 
proceedings). Some other elements whose abundances have been revised due to 
improved atomic data and/or line selection in recent years are Ti, Ni (Reddy et 
al. 2003), Ge (Biemont et al. 1999) and Sr (Barklem & O'Mara 2000). 

Finally, we note that we have not here attempted to homogenise the rather 
diverse collection of photospheric abundances presented in Tabled They have 
not all been derived using the same model atmosphere, departures from LTE 
have only occasionally been accounted for and the methods for estimating the 
attached uncertainties differ between authors. This should be borne in mind by 
the prospective users of this data. 



The solar chemical composition 



11 



2.0 



.0 - 



0.5 - 



0.0 - 



-0.5 







i 



T O 




20 40 60 

Atomic number 



Figure 1. Comparison of photospheric and meteoritic abundances (as mea- 
sured in CI chondrites). The elements which have been analysed using a 3D 
hydrodynamical solar model atmosphere are shown as filled circles while the 
ID-based results with open circles. H, Li and the noble gases all fall outside 
the figure due to significant depletion either in the Sun or in the meteorites, 
which also affects C, N and O 



4. Comparison with the Meteoritic Evidence 

In Table ^ we also present a compilation of meteoritic abundances as measured 
in CI chondrites. The data is largely taken from Lodders (2003) but placed on a 
slightly different absolute abundance scale. It is important to realize that since 
H, like other volatile elements, is highly depleted in meteorites, another reference 
element must be found. The honour is normally given to Si. As a consequence, 
in order to anchor the meteoritic and photospheric absolute abundance scales to 
each other a knowledge of the photospheric Si abundance is required. Since our 
recommended value is 0.03 dex lower than that advocated by Lodders (2003), 
we correspondingly adjust all meteoritic abundances by this amount. 

Fig. ^ shows the difference between the photospheric and meteoritic abun- 
dances. In general the agreement is very good. Some notable differences can be 
easily explained by depletion in the Sun (Li) or in the meteorites (H, C, N, O 
and the noble gases). In addition there are a small number of elements which 
show discordant results: Ga, Rb, Ag, In and W, with CI and Au only barely 
agreeing within the large uncertainties. In the case of In and W one may suspect 
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unidentified blends being the cause of the discrepancies given that the photo- 
spheric abundances are significantly higher than the meteoritic values. For the 
other problematic elements erroneous ^/-values or departures from LTE may be 
to blame. As mentioned above, the new Na abundance does not overlap with the 
meteoritic value. Although the difference is only 0.1 dex, this deserves further in- 
vestigation. The mean difference for the non-depleted elements is 0.01 ±0.06 dex, 
when ignoring the above-mentioned seven elements. 

5. Trouble in Paradise? 

The new solar abundances described here have some important ramifications. 

The first one is evident from the abundances by mass of H, He and metals, 
X, Y and Z, which now become X = 0.7392, Y = 0.2486 and Z = 0.0122 
with ZjX = 0.0165, i.e. much lower metallicity and ZjX than previously 
recommended (e.g. Z/X = 0.0275 from Anders & Grevesse 1989). These are 
the present-day photospheric values with the He abundance taken from Basu & 
Antia (2004). 

The low solar C and O abundances are now in much better agreement with 
those measured in the local interstellar medium (e.g. Andre et al. 2003) as 
well as for nearby B stars (e.g. Sofia & Meyer 2001). Previously the high solar 
abundances apparently already in place 4.5 Gyr ago were inconsistent with the 
low values found in the solar neighborhood today in the framework of Galactic 
chemical evolution. Somewhat surprisingly, the problem appears to have been 
in the solar determination rather than in the analyses of nebulae and hot stars. 

While the new abundances may have solved some long-standing problems, 
it also introduces at least one new, namely for helioseismology. Over the years, 
standard solar interior models which account for element settling have yielded 
an impressive agreement with the observed sound speed in the Sun as measured 
by the solar oscillations. Since C, N, O and Ne are important opacity con- 
tributors, the revised solar abundances for these elements significantly alter the 
interior structure. Indeed the predicted sound speed in these new solar models 
are in much worse agreement with helioseismology (Bahcall et al. 2004; Basu & 
Antia 2004). According to these authors, the only possible culprits are the new 
photospheric abundances or the opacities, having discarded the possibility of 
underestimated element settling in the solar models. The necessary opacity in- 
crease near the bottom of the convection zone is significant, at least 10%. Some 
of this may in fact already have been found (Seaton & Badnell 2004), but more 
than half still remains. Whether the fault lies entirely with the opacity calcu- 
lations or whether our new abundances are in fact underestimated is an open 
question, which urgently needs to be resolved. Experience constantly reminds us 
that one should never claim anything as impossible, but we consider it unlikely 
that our derived abundances are as far from the real values as suggested by 
these helioseismology studies. It would be a remarkable conspiracy of factors if 
this was the case, considering that the atomic and molecular-based abundances 
now finally agree, the lack of any significant trends in derived abundances with 
for example line strength or excitation potential and the almost perfect match 
between observed and predicted line profiles, including their asymmetries. The 
employed lines in our work probe very different parts of the solar atmospheric 
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with greatly different sensitivities to the atmospheric conditions, which implies 
that a great deal of fine-tuning would be necessary to simultaneously bring the 
C, N and O abundances up by some 0.2 dex, if at all possible. As mentioned 
above, the new abundances also removes the special nature the Sun otherwise 
would have in comparison with its neighborhood. 

6. Concluding Remarks 

Although progress has been made in refining our knowledge of the solar chemical 
composition, much work still remains. An obvious way forward is to extend 
the 3D analysis presented here to the remaining elements in the periodic table. 
Furthermore, the lack of detailed non-LTE line formation calculations is a major 
cause of concern, not only for the solar abundance determinations but also for 
much of the stellar work. Equally important, however, is continued investment 
in improving the necessary atomic and molecular data. Although the impact 
may not be as profound as the recent surprisingly large revisions of the solar 
metal content, efforts along these lines would nevertheless be worthwhile. 
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